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ABSTRACT 

We report on average sub-surface properties of pre-emerging active regions as compared to areas 
where no active region emergence was detected. Helioseismic holography is apphed to samples of 
the two populations (pre-emergence and without emergence), each sample having over 100 members, 
which were selected to minimize systematic bias, as described in Paper I (Leka et al., 2012). We 
find that there are statistically significant signatures (i.e., difference in the means of more than a 
few standard errors) in the average subsurface fiows and the apparent wave speed that precede the 
formation of an active region. The measurements here rule out spatially extended flows of more than 
about 15 ms~^ in the top 20 Mm below the photosphere over the course of the day preceding the 
start of visible emergence. These measurements place strong constraints on models of active region 
formation. 

Subject headings: Sun: helioseismology - Sun: interior - Sun: magnetic flelds - Sun: oscillations 



1. INTRODUCTION 



As discussed in Leka et al. (2012) (Paper I), the mech- 
anisms behind the formation of solar active regions ( AR) 
are not known. Possibilities include magnetic flux tubes 
rising essentially intact from the base of the convection 
zone (for a review see Fan 2009). Alternatively, near- 
surf ace effects could dom mate the formation mechanism 
(e.g. Brandenburg|200"5 ) ; hybrid scenarios are also a pos- 



sibility. More generally, understanding of the formation 
of active regions may lead to a better understanding of 
the solar dynamo. 



Local h elioseismology (Gizon & Birch 2005 Gizon 



et al. 2010, ) is among the tools that potentially could 



be used to determine the subsurface dynamics associ 
ated with active region formation. While helioseismic 
analysis samples the state of the plasma below the visi- 
ble surface, the interpretation of the results is not always 
straightforward - especially when the region sampled is 
quickly changing in context, such as an emerging mag- 
netic flux region. 

Numerous studies have attempted to detect - and thus 
characterize - the signature of one or two active regions 
at their earliest st ages o f formation. In one of the e arli- 
est studie s, Braun|(1995) applied Hankel analysis (Braun 
2) to South Pole observations of the lor- 
A AR5247. In the few days before the 
region's sunspot was reported, the Hankel analysis de- 
tected negative phase shifts (i.e., reduced wave speed). 
The cause of the phase shifts is not known, and the sup- 
p orting magnetic me asurements were very s parse. 



[ Chang et al. ( 1999[) used acoustic imaging ( Chang et al. 
1997') applied tVMDI 



and also TON data ('(Jhou et al. 



1995) to study the emergence of NOAA AR7978. Based 
on an analysis of the focus-depth dependence of observa- 
tions with a two-day resolution, they suggested that they 
had detected a magnetic flux concentration moving up- 
wards towards the photosphere. However, the detected 
positive phase shifts that are co-spatial with the emerg- 
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ing active region are also co-temporal with the appear- 
ance of su rface fl ux. 

Jensen et al. ( 2QQT|) us ed time-distance helioseismol- 
ogy ([Puvall et ar"T593) to study two active regions 



that formed on 11 January 1998 (NOAA AR8131 and 
AR8132) and found changes in the subsurface wave 
speed that developed at about the same time as the mag- 
netic signatures of the regions were seen in MDI magne- 
tograms. There was no apparent change in subsurface 
st ructure before the su rface magnetic features were seen. 



Hartlep et al. (2011 ) measured the acoustic power from 
MDI observations before and during the emergence of 
AR 10488. This study suggested a decrease in acoustic 
power in the 3 to 4 mHz band and a reduction in subsur- 
face wavespeed before signiflcant magnetic flux is seen at 
t he photosphe re. 



Kosovichev 



(2009) used time-distance helioseismology 
to study two emerging active regions (NOAA AR8167 
and AR 10488). These case studies inferred increases 
in the subsurface wave speed associated with flux emer- 
gence, however there was no clear indication that these 
increases preceded the appearance of magnetic fl ux at 
the photosphere. |Zharkov & Thompson (2008) also 
used time-distance helioseismology in a case-study of the 
emergence of NOAA AR 10790; sound-speed increases 
were also found with (but not prior to) the appearance 
of surface fleld. 

f 

th 



Komm et al. 
to infer 



[iSHSf 



(|2008| used ring-diagram analysis ( |Hill 
le subsurface flows; it was applied to 



a lew active regions undergoing flux emergence: newly- 
emerging regions (including AR 10488) were compared to 
older but resurging active regions. The analysis, which 
included a comparison with magnetic fleld in the target 
areas, found indications of upflows prior to emergence 
that were followed by downflows. 

In an extension of the ab ove study to one of the few 
truly statistical approaches, |Komm et al.| ( [20091 [20TT| 
then used the same ring-diagram analysis approach on 
a very large sample of active regions characterized by 
the evolution of the magnetic fleld (from increasing to 
decreasing) . Key to this study was the "control" dataset 
of as many "quiet" areas. This study found that growing 
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active regions show a preference for upflows of a fraction 
of ms~^, with this effect most apparent at the depth of 
about 8-10 Mm and reversed below this depth. 

There are indications that the details of particular data 
analysis methods may be very important in determin- 
ing the helios eismic signal s seen for emerging active re- 
gions. Ilonidis et al. (2011) showed several cases (includ- 
ing AR 10488) where a travel-time reduction of order 10 s 
for measurements sensitive to a dep th 60 Mm wa s seen to 
precede the emergence of an AR. Braun (2012) showed 
that this result is not reproduced (on the same active 
regions) using helioseismic holography. 

It is difficult to get a consensus picture of the results, 
likely for a variety of reasons. While most studies re- 
port some detection of a signal as sociated with t he ap - 



pearance of an active region, only Komm et al. 
Hartlep e ^aTI (|2011|), |Komm et aT.| ( |2009 



and 



( 20081), 



lomm 

7T 



et al. ( 20lT| ) include aiiy comparison to a ''contro 
pie by whicn to evaluate the detections, although in most 
cases the criteria for this evaluation are not discussed in 
detail. The presence of surface fields may complicate the 
interpretation of any 'pre-emergence' signal. 

While mention of the surface field evolution was in- 
cluded in most studies, a detailed accounting of the pre- 
sense (or absence) of surface-field for the full time period 
used for the helioseismology data analysis was rarely in- 
cluded. The studies cited above are predominantly case 
studies, with the analysis performed at (or over) varying 
times in their target regions' evolution. Case studies are 
enlightening and critical f or guidance - but as discussed 



Paper I and shown in Ilonidis et al. (2011), regions 



emerge with different rates and into ditterent surface con- 
texts; this points to the need for statistical studies, but 
also reveals why few strong results have emerged f rom th e 
large-sample studies to date (Komm et al. 2009] |20lT ). 

Numerical simulations have been used to model the 
emergence of magnetic flux t hrough the last few tens of 



Mm below the photosphere. Stein et al. (2011) used a 
radiative magnetoconvection simulation in which hori- 
zontal magnetic field was advected into the simulation 
domain by upflows through the lower boundary (20 Mm 
below the photosphere). For the case of 10 kG field 
strength, the magnetic field emerges through the photo- 
sphere in about 32 hours. In this case, the field has only 
a weak effect on the convection and the emergenc e pro- 
cess is largely controlled by the convection. Cheung et al.| 
(^2010), also using a radiative magnetoconvection simu- 
lation, studied the emergence of a semi-torus of twisted 
field. In this case, the magnetic structure was forced 
through the lower boundary with a vertical velocity of 
1 km s~^ and took about 2 hours to rise to the photo- 
sphere from a depth of 7.5 Mm. In these simulations, 
vertical flows of about 0.5 km s~^ and diverging horizon- 
tal flows of order a few km s~^ are associated with the 
flux emergence at the photosphere. These two simula- 
tions show scenarios for flux emergence with very differ- 
ent observational implications. 

In the present study, we combine a statistically- 
significant sample and a control group , with a detailed 
analysis of the surface magnetic field ( |Paper ip , and ex- 
amine the differences in pre-emergence behavior in an 
average sense. We undertake a more detailed statistical 
analys is of the two popula tions in the next paper in this 
series ( Barnes et al.||2012 , Paper III). It should be noted 



that the analysis used herein (see q3| relies on "raw" 
travel-time shifts rather than inversions used by most 
of the studies cited above; as such, we avoid potential 
complications and uncertainties of the inversions, while 
sacrificing details about the depth of any perturbations 
detected. Nonetheless, we show that on average there 
are statistically significant sub-surface flows and changes 
in wave speed that precede AR emergence; these results 
place strong constraints on models of active region for- 
mation, even when individual case studies show no clear 
signature of emergence. 

The layout of the paper is as follows. In ^ we review 
the data sets that are used in the helioseismic measure- 
ments. In ^ we detail the application of helioseismic 
holography to these data sets. A statistical summary of 
the results is shown in We discuss the main results 
and outline some possibilities for future work in ^ 

2. THE DATA 

The GONG Doppler data cubes for active regions 
before their emergence (we will refer to these as pre- 
emergence, or "PE" cases) and for quiet- Sun control 
cases that are not associated with AR emergence (non- 
emergence, or "NE", cases), accompanied by the as soci- 
ated M DI magnetograms are described in detail by |Pa-| 
per I[ To briefly review the procedure described there: 
we selected a sample of 107 pre-emergence cases and 
the same number of non-emergence cases. These sam- 
ples of the two populations (PE and NE) were selected 
to have well-matched distributions in disk location (lat- 
itude/CMD) and time (within the solar cycle), to avoid 
biases in the seismology due to projection and instru- 
mental effects. For the PE cases, a refined emergence 
time to was defined as when the change in the total ab- 
solute flux (as measured from MDI) reached 10% of the 
maximum increase detected within a 3-day window of its 
nominal (NOAA-defined) emergence time. The GONG 
data cubes were 1664 minutes long (one "GONG-day"), 
spanning from 1648 minutes before to to 16 minutes af- 
ter. These 1664 minute cubes were divided into five time 
intervals of 384 minutes each, with an overlap of 64 min- 
utes between time intervals. The time between the start 
of each time interval is 320 minutes. The cadence of the 
GONG Doppler data is one minute per image. 

As described by Pap er Ij each o f the Doppler images 
is Postel projected (Pearson 1990) on to a map with a 
scale of 1.5184 Mm pixel" ^ and a size of 256x256 pixels^. 
Figure [l] shows the mean power spectrum of the NE re- 
gions. The ridges in the power spectrum are visible up 
to about /c = 1.5 radMm"^. 

Potential field extrapolations from MDI magnetograms 
are used to estimate the r adial com ponent of the mag- 
netic field (for details see Paper I). These estimates of 
the radial magnetic field are then projected into the same 
map coordinates as the GONG Dopplergrams and aver- 
aged over the same time intervals. The result of this 
procedure is a set of maps of radial magnetic field esti- 
mates, one map for each time interval for each PE and 
NE case. 

3. HELIOSEISMIC HOLOGRAPHY 



Helioseismic holography (|Lindsey & Braun 2000) is a 
tool that uses measurements of solar oscillations to infer 
subsurface conditions; it is very similar to time-distance 
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Figure 1. Mean power spectrum of the GONG Doppler data 
cubes for all of the NE regions. The ridges are visible up to about 
horizontal wavenumber k = 1.5 rad Mm"-*^ . The central phase 
speeds for the filters used in this work (we use the filters from 
[Couvidat et al.|2005| > are shown as dashed lines (filter TDl has the 
lowest phase speed and TDll has the highest). The width of the 
filters is similar to the distance between filters. 



filter 


inner radius ( Mm) 


outer radius (Mm) 


depth (Mm) 


TDl 


3.7 


8.7 


1.4 


TD2 


6.2 


11.2 


2.2 


TD3 


8.7 


14.5 


3.2 


TD4 


14.5 


19.4 


6.2 


TD5 


19.4 


29.3 


9.5 


TD6 


26.0 


35.1 


11.4 


TD7 


31.8 


41.7 


13.3 


TD8 


38.4 


47.5 


15.7 


TD9 


44.2 


54.1 


18.2 


TDIO 


50.8 


59.9 


20.9 


TDll 


56.6 


66.7 


23.3 



Table 1 

Table of annulus size and lower turning point depth, taken from 

Table 1 of Couvidat et al. (2005). The first column shows the 
filter name. I'he remaining three columns show the inner annulus 
radius (Mm), outer annulus radius (Mm), and the lower turning 
point depth (computed at the central phase speed of each filter). 



helioseismology ( Duvall et al.||1993 ). See Gizon & Birch 
d2005) for a detaileT review ot helioseismic holography. 
Here, we apphed surface- focusing hehoseismic hologra- 
phy to each of the five time intervals of each of the PE 
and NE data cubes. 

The basic steps in the analysis are: 1) track and Postel 
project GONG Dopplergrams (described in the previous 
section), 2) apply phase-speed filters, 3) compute local- 
control correlations , and 4) measure travel-time shifts. 
Phase-speed filters (Duvall et al. 1997| ) isolate waves with 
particular ranges in lower turning points. We used filters 
1 through 11 (here denoted as filters TDl to TDll) as 



described in Table 1 from Couvidat et al. ( |2QQ5 ). These 
filters cover the range in sampling depth from about 
1.4 Mm (filter TDl) to about 23.3 Mm (filter TDll), as 
shown in Table [T] Each phase-speed filter leads to a sepa- 
rate filtered data cube (i.e., filtered time series of Dopp- 
lergrams). From each of these filtered data cubes, we 



then computed local-control correlations (e.g., Lindsey 



Braun||2QQQ[ these are analogous to time-distance cor 
relations) for both center- annulus and center- quadrant 



geometries (see e.g., Duvall et al.||1997 , for a description 
of these geometries), 'i'he annulus sizes are shown in Ta- 
ble [l] and were taken from Table 1 of ICouvidat et ah] 
( 2005 ) . Finally, we measured travel-time shifts from the 
center- annulus and center- quad rant local- control correla - 
tions using the phase method of ] 
We emphasize that due to the 



Lindsey fc Braun| ( |2000[ ). 
racking procedure, most 



of the effect of solar rotation is removed from the mea- 
surements shown here. In addition, as described in Q 
we remove smooth fits to all of the maps of travel-time 
maps. As a result, the travel-time maps shown here do 
not include large-scale effects (e.g., differential rotation 
and meridional fiow). 

Throughout this paper, we use x and y to denote the 
coordinates in the Postel projection geometry, with x 
increasing in the direction of rotation and y increasing 
to the north. We use 5rx to denote east-west travel- 
time differences (i.e., the difference in travel times be- 
tween east-going and west-going waves) and Sry to de- 
note north-south travel-time differences. These travel- 
time differences are mostly sensitive to east-west (north- 
south) fiows, with a negative 5rx {^Ty) indicative of a 
fiow to the west (north). 

Following the usual convention, we denote the outgoing 
(i.e., center-to-annulus) travel-time shift (shift relative to 
the average) by Sr out and the in-going (i.e., annulus-to- 
center) travel-time shift by Sr^^. From these we construct 
the "out minus in" travel-time shift = 5r o^t — ^^in 
and the mean travel-time shift Sr mn = [^^^out + ^^^in] /2. 
The one-way center- annulus travel times Sri^ and ^Tout 
are sensitive to the isotropic wave speed and also to 
vertical fiows and converging/diverging horizontal flows. 
These two effects (i.e., change in isotropic waves speed 
and the presence of flows) are approximately separated 
by ^Toi, which is sensitive mostly to diverging / converg- 
ing horizontal flows and vertical flow, and Srmn which is 
more sensitive to changes in the isotropic wave speed. 
The sign conventions are such that Srmn < is in- 
terpreted as a signature of increased wave speed, and 
5roi < is the signature of a diverging flow. 

From the quantities 5rx and 5ry we computed prox- 
ies for the vertical component of the flow vorticity (de- 
noted vort) vor = dxSry — dySrx and for the horizontal 
flow divergence (denoted div) as div = dxSrx + 9y6Ty. 
The derivatives were computed in the Fourier (horizontal 
wave vector) domain. 

4. RESULTS 

The analysis is presented in three stages: flrst, indi- 
vidual examples; second, averages over all samples (sep- 
arately for the two populations) but retaining the spatial 
information; and third, an analysis of the spatial average 
over a small central region, and its temporal evolution, 
again contrasting the averages over the two populations. 
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4.1. Individual Active Region Examples 

Figure [2] shows, for the case of the emergence of 
AR 9729, the Sr ^nn maps for filter TD5 (lower turning 
point of about 9.5 Mm) together with estimates of the 
radial magnetic field obt ained fro m potential field ex- 
trapolations (details are in [Paper "l| from MDI 96-minute 
magnetograms. In this case, the emergence takes place 
in a region with strong nearby magnetic fields. These 
pre-existing surface magnetic fields have corresponding 
negative features (i.e. increased apparent wave speed) in 
the maps of Srmn- Using the NE cases, we estimate that 
the noise level in the smoothed maps (for maps with duty 
cycle of better than 70%) of Srmn for filter TD5 is about 
0.9 s. Therefore, many of the weaker features seen in 
the travel-time shifts in Figure |2] are likely due to noise. 
With this noise level in mind, there is no apparent evo- 
lution with time. Notice that this figure covers roughly 
the 24 hours before emergence, and as a result the de- 
velopment of the active region is not seen in either the 
magnetograms or the travel-time maps. This is a typical 
example of the PE cases. 

Figure |3] shows another example, this time for the pre- 
emergence time evolution of NO A A AR 10488, that was 
also studied by Komm et al. ( 2008|), Kosovichey fc Du- 



i^('20Q8^,'KQSQv khey^^ et al.|( |:^011p iSa' 

rionidis et al. (20lT|. in our analysis, there is no discern- 
able signal greater than ±2 s in the first four time inter- 
vals (covering approximately —27 hr < to < —5 hr) in 
the raw or smoothed Srmn maps; the raw and smoothed 
magnetic flux maps are similarly bare. In fact, there are 
fewer areas of nearby strong field and spatially-associated 
Srmn signal than in the previous example (Figure |2|. 

In this case, a feature with ^ 50 G magnitude appears 
in the average magnetogram for the last time interval. 
This is not strictly unusual; the automated definition of 
emergence time that we have used, while refined earlier 
by hours or days from the NOAA-assigned emergence 
time, may still allow magnetic field at the surface prior 
to the "emergence ti me" to o f 11:11 UT on 26 October 
2003 in this case (see Paper I for details). Here we note, 
however, that less objective definitions of emergence lead 
to less repeatable results; flux emergence does not al- 
ways follow a standard template for temporal evolution. 
Nonetheless, the signal is weak in this last time interval 
(it is an average over ^ 6 hr) and has no corresponding 
features in the maps of Sr mn^ though this may be due 
mostly to noise (which we expect, based on the statis- 
tics of the NE cases, to be about 0.9 s in the smoothed 
travel-time maps for the filter shown in Fig. [3|. 

Of the other studies of this region, all use a different 
emergence time; here, we summariz e the results relative 
to our definition of emergence time. Komm et al. (2008) 
begin their analysis on 27 October z003, and thus there 
is no possibility for a direct comparison of the results, 
as this is the day following our analysis. [Kosovichev &[ 
Duvah (2008); Kosovichev ("2009") first see wave speed 



perturbations in data from an 8 hr time interval centered 
at approximately one hour after to, and describe this as 
a pre-emergence signature with growth of the magnetic 
fiux starting at approximately to + 9 hr. Thus, the results 
are similar to this study in that no clear subsurface signal 
appears when only data from at least a few hours prior 
to to are used in the analysis. 



Hartlep et al. (2011) state that magnetic field starts to 
appear at the surface approximately 2 hr before to, while 
significant magnetic fiux appears at around to- Their 
finding of a change in the acoustic power in the 3-4 mHz 
frequency range, averaged over 128 minutes, centered ap- 
proximately two hours before to is likely related to the 
weak magnetic field seen in the last time interval of this 
study. 



Finally, Ilonidis et al. (2011) report a reduction in the 
mean travel-time of about 16 s at a depth of about 60 Mm 
(much deeper than the 9.5 Mm shown here) from an 8 hr 
data set centered approximately 7.5 hr before to- The 
difference in the depths conside red aga i n mak es a direct 
comparison impossible, but see Braun (2012) for a dis- 
cussion of using helioseismic holography to examine the 
same depth range. 

Even for this most frequently studied active region 
emergence, little direct comparison is possible due to dif- 
ferent time intervals and depth range s covered. However, 
in the mo st c omparable s t udies of [Kosovichev fc DiT 
vail (2008) and Kosovichev (2009), the results are similar 
to the results presented here when the difference in the 
stated time of emergence is accounted for. 

4.2. Averages for the Two Populations 

We now present the results of averages over all mem- 
bers of a sample (e.g., the average ove r all NE cases) 
with duty cycle greater than 70% (c.f.. Paper I Table 
3; all averages in this paper will use this constraint on 
the duty cycle). Averages taken in this manner are in- 
dicated by angle brackets, (•). For the travel-time maps, 
large-scale spatial variations resulting from the Postel 
projection have been removed prior to averaging, using 
a second order polynomial fit in the x and y map coor- 
dinates (the functional form of the fitting polynomial is 
f{x^y) = aix'^ + a2y'^ + ci?>xy + a/^x + a^y + a^). For 
this analysis, all maps were smoothed with a Gaussian 
of FWHM of 9.99 pixels. In the case of the PE sam- 
ples, the emergence site was not initially at the center of 
every data cube. To improve the signal/noise ratio for 
this analysis, we aligned all of the travel-time map s (and 
magnetograms) before averaging. As described by Paper 
I, the alignment was done by first computing, for each 
PE case, the centroid of the pixels where the change in 
time averaged between time interval (about 27 to 
21 hours before emergence) and the average of Br com- 
puted over the 4 to 10 hours after emergence was more 
than 30% of the maximum. 

Figure [4] shows the time-evolution of {Srmn)^ (Stx), and 
(Sry) maps over all of the NE cases, for filter TD5. Also 
shown is the evolution of the corresponding average of 
the unsigned radial magnetic field. The average magnetic 
field is overall weak and unchanging over these averages. 
There are regions where the average magnetic field is 
above 20 Gauss, generally near the edges of the map. 
Nearby plage field was pern iitted in the definition of a 
"non-emerging" region (see Paper I[ for a discussion). 
The maps of (Srmn)^ {Stx)^ and (^r^y), all for filter TD5, 
show no features stronger than ±ls, and little spatial or 
temporal coherence between the time intervals. 

Figure [5] shows the time-evolution of the same aver- 
aged quantities for the PE samples (spatially co-aligned 
as described above). The difference between these re- 
sults and those shown in Figure [4] is striking. There is a 
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Figure 2. Example, for NOAA AR9729, of the time evolution of the (5rmn and of the surface magnetic field over the 24 hr prior to 
the appearance of the region's emerging magnetic fiux. Time advances from left to right in steps of 320 minutes; the mid-interval time 
relative to to is labeled at the top. Top row: maps of (5rmn for filter TD5 (corresponding to a depth of about 9.5 Mm) for each of the 
time intervals. The last column covers — 368min < to < 16min, little if any active-region fiux emergence should be visible. Second row: 
same as the top row after smoothing with a Gaussian with full- width at half maximum (FWHM) of 9.99 pixels (about 15 Mm). Notice the 
change in color scale; the spatial smoothing reduces the amplitude of the fiuctuations in (5rmn- Third row: the evolution of the unsigned 
radial magnetic field strength. Fourth row: same, but after the same smoothing as was applied to the maps of (5rmn; again, the color 
scale changes because of the reduced amplitude of the fiuctuations. Notice that negative trav el-ti me shifts appear co-spatial with the areas 
of strongest magnetic field. The black circle in all panels shows the averaging region (see § |4.3| ) surrounding the area where the AR will 
emerge. The inset (top left panel) indicates the size of the annulus used in the holography measurements for filte r TD 5. Notice the higher 
noise level in 5rmn maps for the second time interval; this is caused by a low (60%) duty cycle. As discussed in ^4.2[ such low duty cycle 
time intervals are excluded in later analysis. 



persistent feature of order -1 s in the maps of ((^Tmn) for 
filter TD5, with the strength of the signal increasing as 
to approaches. This signal appears to be associated with 
a feature in the average magnetogram in the same loca- 
tion. In general, we see correspondence between strong 
features in the average radial magnetic field strength and 
features in the mean travel-time shifts. One possibility is 
that the surface magnetic field is the direct cause of the 
travel-tiine shi ft (i.e., the showerglass effect of [Lindseyj 
& Braun||2005|). We undertake a statistical study of this 



beginning a day prior to to- For the first four time inter- 
vals, there are hints of anti-symmetric features in these 
maps; a fiow converging towards the emergence site at 
approximately 15 ms~^ might give such a distribution 
of travel times. The maps for the last time interval are 
more complicated; in particular the map of {^Tx) shows 
what might possibly be the signal of t he magnetic region 
moving in the prograde direction (e.g., Zhao et al.|2QQ4 ). 



possibility in the next paper in this series (Barnes et al. 
|2Q12 , Paper III). 

We also see (Fig. |5| that the average radial field 
strength at the emergence site also increases as to ap- 
proaches. While the definition of to allows nearby fiux 
to be present prior to to, and may be uncertain to a few 
hours (see Pa per I) , a signal can be seen at least a day 
before emergence. There is a clear preference for sur- 
face field to be located at th e emerge nce site prior to 

— 



and 



Paper III 



significant fiux emergence (see [Paper I 
further discussion) 

Figure [5] also shows there are persistent features in the 
maps of (OTa^) and ((^r^y) with amplitudes of about 2 s 



4.3. Spatial Averages for the Two Populations 

Motivated by Figure [5j and to focus more closely on 
the site of fiux emergence, we next present spatial aver- 
ages. The averages are computed, for each sample in the 
PE and NE populations, over disks of radius 45.5 Mm 
centered on the emergence location (as de fined usin g the 
centroid mentioned above, c.f. details in Paper I), and 
indicated in Figures |2]-[5j We denote these initial spatial 
averages with an over-bar, e.g., Stx- As discussed ear- 
lier, data cubes with duty c ycle less than 70% are omit- 
ted from this analysis (c.f.. Paper I Table 3). We then 
compute averages over the entire PE and NE samples, 
separately for each time interval, and the standard error 
in the mean (cr/y^, where a is the standard deviation 



6 



60 Mm 



60 Mm 



-19.2 hr 



-13.9 hr 



-8.5 hr 



-3.2 hr 







I 







-5 ^ 



-10 
4 

2 ^ 

GO 



_2 ^ 



1600 



400 § 

bJO 



200 



I"0 
250 

200 ^ 

150 i 

bJO 



100 . 

150 




Figure 3. Same as Figure |2] but for AR 10488. As in Figure [2] notice 
for the third and fourth rows. In this case, there is no clear signature 
before emergence. 

and n is the number of samples included). The results, 
described further below, are depicted in a multi-panel 
presentation that includes the parameters, all filters, and 
the temporal evolution of both PE and NE averages and 
their errors within each filter/parameter combination. 

Figure [g] shows the time evolution of (^Tmn), {^^x)-, and 
{5Ty). The figure also shows {Stx cos 0) and (Jr^ySin^); 
these quantities are the sample average (angle brackets) 
of the spatial averages (over-bar) of either Stx weighted 
by cos^ or Sry by sin^, where is the angle measured 
counter-clockwise from the direction of solar rotation 
(the +x direction). 

Focusing first on the measurements ((5rmn), there ap- 
pear to be clear and sustained differences between the PE 
and NE cases for ah but the shallowest filters (TDl, TD2, 
and TD3), and the deepest filter (TDll). The sense of 
the difference is generally that {6rmn) is more negative 
for the PE cases than the NE cases. The NE cases are 
very temporally consistent between filters, and in most 
filters the difference between the NE and PE samples in- 
creases as time approaches the emergence time to, with 
the PE population results increasing in (negative) mag- 
nitude. The amplitude of this effect is of order tenths 
of a second. This may perhaps be a consequence of the 
surface magnetic field and is discussed in detail by [Paper] 

InD 

Regarding the fiow measurements (Stx) and (Sry), 
there is no clear difference between the NE and PE pop- 
ulations. The anti-symmetric averages {5rx cos 0) and 



that the color scale is different for the first and seconds rows, and 
of emergence in the magnetic field or 6rmn maps in the 24 hours 



{dry sin 0)^ however, show substantial differences (e.g. of 
between two and three standard errors for filter TDS) for 
filters TD2-TD5 (depths of 2.2 Mm to 9.5 Mm). These 
differences are due to the anti-symmetric features seen in 
Figure [s] Notice that in all of these filters (TD2-TD5) 
these features are seen even at 24 hours before emergence. 

Figurej7[ shows the time evolution of ((^rin), (^^Tout), 
(^''"oi), (div), and (vor). The parameters (Srin) and 
(^'^out) show fairly consistent differences between the PE 
and NE samples. These appear to be roughly consistent 



with the change in ((5rmn) seen in Figure [6j Similarly, 
there is a trend that (Sroi) is systematically larger in the 
PE than the NE cases, again this is qualitatively consis- 
tent with the maps of (Sr^) and (Sry) from Figure [s] 

Other results are less pronounced. For example, in 
(div) there is a weak difference between the NE and PE 
samples for filters TD4, TD5, and TD6; the sign is con- 
sistant with a converging fiow. Other PE/NE differences 
are either fieeting or at the one-<j level. For example, the 
average divergence (div) shows a strong positive feature 
for filter TDl in the second time interval, this feature is 
not persistent in time and may be noise. In the first time 
interval, there are differences, at the one-a level, between 
the (vor) for the PE and NE cases for filters TD2, TD3, 
TD5, and TDll. There are also one-a differences for the 
last time interval for filter TD2, the second time interval 
for TD7, and the last time interval for TD9. 

In the interest of recognizing spurious results, we note 
the following. Assuming that the variances of two dis- 
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Figure 4. Evolution (left to right, times labeled as Figure[2| of {\Br\), {Srmn) ,{^Tx) , and (Sry) maps for filter TD5, of the average for each 
time interval taken over all NE samples with duty cycle great than 70%. No clear and persistent features are visible, and no features are 
present that are distinctly greater than the noise. The black circle in each panel indicates the area that is spatially averaged in subsequent 
analysis. 



tributions are equal (e.g., PE and NE; the assump- 
tion is reasonable if the variance is predominantly due 
to noise), simply by chance alone approximately 16% 
of the results should have the difference between the 
means greater than the sum of the standard errors (i.e., 
\xi — X2\ > (Ti/s/ni -\- cr2/y^, where a is the standard 
deviation, and as mentioned above, a j ^Jn is the stan- 
dard error shown in Figures [g] and [t]) . This translates 
to roughly nine 1 — a results per parameter that are ex- 
pected to be purely spurious. For the difference being 
twice the sum of the standard errors, the expectation 
shrinks to 0.5%, or less than one result per parameter, 
and for the difference being three times the sum, the ex- 
pectation would be a nearly negligible 3 x 10~^%. For 
the case of (vor), there are seven such points. Hence, it 
seems possible that the differences we have seen in (vor) 
and (div) as well as some of the other differences between 
the PE/NE samples are simply the result of noise. Con- 
versely, the many large differences in ((^rin), for example, 
are very likely real. 

4.4. Results for the Ultra- Clean^' Pre- Emergence 
Sample 



In [Paper! we describe a subset of the PE sample which 
we deemed ultra-clean" . These 11 hand-selected emerg- 
ing active regions all eventually reached a size of at least 
70 /iH and displayed a monotonically increasing flux his- 
tory after a distinct and unambiguous emergence time 
to, appearing into an otherwise weak field area. 



Figure [s] shows (|5r|), (^Tmn), i^'^x)^ and (Sry) aver- 
aged over this subsample of the PE cases. As in the 
average over all PE cases, the strongest magnetic feature 
appears at the emergence site in the last time interval. 
Also as in the average of all PE cases, there is an appar- 
ently associated feature in the average of 6rmn- Inter- 
estingly, a similar feature in the -13.9 hr interval has no 
corresponding feature in (l^rl), and there is no (^Tmn) 
feature in the interim -8.5 hr interval. Due to the smaller 
number of regions, the noise is higher in these averages 
(we expect the noise to be about 0.9/\/l0 ~ 0.3 s in the 
smoothed maps of {Srmn))^ thus, it is not clear if lack 
of features in the average Srmn map at the earlier times 
is significant. The {drx) and (Sry) maps are also noisier 
than the corresponding maps in Figure [5j as expected. 
There is still a weak suggestion in some of the time in- 
tervals of anti-symmetric features in the average maps of 
Stx and Sry (again as in Fig. l5|. Still, it is worth em- 
phasizing that these 11 hand-picked regions presumably 
represent the best-case (i.e. relatively simple emergence 
into quiet Sun) scenario for detecting helioseismic signals 
in the day prior to an active region's emergence. 

Figure j9]_shows the spatial and ensemble averages 
(Jrmn), {^'Tx)^ {^'^y)i (^'?"ic cos^), {dTy^iiiO) for the sub- 
sample of eleven PE regions, along with the same average 
over the NE regions shown in Figure (6] The results are 
less clear than what is presented in Figure [6) due to the 
much smaller sample size for the PE population (about 
nine instead of 90, after accounting for time-intervals 
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Figure 5. Same as Figure [4] showing (top: bottom) {\Br\)^ (<^Tmn}, {^^x), and (Sry) but for the average over all of the PE samples 
(with duty cycle great than 7D%) again for filter TD5. The features in the maps of {Srmn) are spatially correlated with those in the maps 
of {\Br\). For the first four time intervals, there are anti-symmetric features in the travel-time differences (Stx) and (Sry); these features 
might be hints of a converging fiow. The map of (Srcr) for the la st time intervals shows what might be evidence of a weak prograde fiow 
associated with the surface magnetic field (e.g. [Zhao et al.|2004| ). These travel-time differences are of order two seconds, which for filter 
TD5 corresponds to a fiow of about 15 m s"-*^. 





with low duty cycle, see Table 3 of Paper I). For this 
subsample of pre-emergence active regions, no difference 
between the (drmn) is seen. This may be because in the 
hand-selected subsample of PE cases, there is very little 
nearby surface magnetic field except in the last time in- 
terval (see Fig. |8|. There is no apparent signal in (Stx) 
or (Sry). The differences between the PE and NE popu- 
lations can, however, still be seen in the anti-symmetric 
averages {Stx cos 6), {Sry sin 6) for filters TD2-TD5. It is 
noteworthy that these differences remain, even when the 
differences in {6Tmn) are not visible. This suggests that 
there are different physical mechanisms responsible for 
these two effects. 

Figur e \iO\ shows the spatial and ensemble averages 
(^^Tin), ((^Tout), (^'^oi), (div), and (vor) for the subset of 
eleven PE regions, following Figure [7| For the variables 
(^'Tin), (^Tout)^ (^'^oi), and (div) there are no statistically 
significant differences in the means for the two popula- 
tions. For the case of (vor), there is still a trend to see 
differences at up to the two-a level for the first time in- 
terval, with the sign of the effect depending on the filter. 

For these results, we note again that the standard er- 
rors for the subset of the PE sample (Figuresj9]and 10) 
are much larger than those shown for the full PE sample 
in Fig.[7| In addition, due to the smaller sample size, the 
chances of spurious differences at and beyond the 1 — a 
level are significantly larger at approximately 21% of the 



results. This translates to roughly twelve 1 — a results 
per parameter that are expected to be purely spurious. 

5. DISCUSSION 

We have applied helioseismic holography to 107 pre- 
emergence active regions and an equal number of areas 
where no active region emergence occurred. The sam- 



ple of emerging active regions, as described in Paper I 
spans a wide range of eventual sizes and temporal evo- 
lution profiles of Br] in addition, these active regions 
emerge into a variety of magnetic contexts. We present 
single-region examples, but concentrate herein primar- 
ily on averages over the samples in the two populations, 
and on spatial (and ensemble) averages over the central 
emergence area. 

We have found differences in the average seismic signa- 
tures between these two populations. The pre-emergence 
regions show a ^Tmn that is reduced by several tenths 
of a second (in all but the shallowest layers), compared 
to the non-emergence population. There are persistent 
anti-symmetric features with magnitudes of up to 2 s in 
Stx and dry in the shallow filters. These features may 
qualitatively suggest a converging flow of order 15 ms~^. 
No clear statistically significant differences between the 
average pre-emergence sample and no-emergence sample 
were found in the vorticity and only hints seen in the 
divergence. 

However, two caveats must be included in the above 
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Figure 8. Same as Figure 14] but for an average over the eleven hand-selected "ultra-clean" PE regions. In this case, the (|^r|) < 20 G at 
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in Figurels] as the average is over a much smaller sample of regions (roughly nine instead of ninety, see Table 3 of Paper I). There is a very 
weak signal in {Srmn) at the time intervals centered at -13.9 hr and -3.2 hr. 



summary. First, there exists a weak but persistent av- 
erage magnetic signal at the emergence site prior to de- 
tectable flux emergence; it strengthens as the emergence 
start approaches, as does the magnitude of the co-spatial 
(5rmn deficit. The causal relationship between these fea- 
tures is not clear: does the change in wave speed and the 
apparent fiow result from the surface magnetic field? has 
the magnetic field collected at the emergence site because 
of the converging fiows? 

Another possibility is that this is another example of 
a bias in the selection of the samples used in this in- 
vestigation, and our results are a consequence of emer- 
gence happening preferentially at the boundary between 
supergranules. Our non-emergence sample has no pref- 
erential location compared to supergranules, but if the 
emergence locations are preferentially centered on the 
boundary between supergranules, this would result in 
fiows towards the emergence location. The fiows would 
not then be converging on the emergence location, but 
rather the emergence location would be between neigh- 
boring diverging fiows from the supergranules. 

The depth of the antisymmetric features in (Stx) and 
(Sry) approximately corresponds to the typical depths 
to which supergranular fiow is detected, ^ 5 Mm, and a 
typical supergranule lifetime of 1-2 days is at least as lon g 
as the time considered here (Rieutord & Rincon 2010|. 
Further, the presence of surface magnetic field prior to 
emergence could also be a result of the tendency for mag- 
netic fiux to concentrate in the boundaries between su- 



pergranules. Finally, this might also account for the com- 
paratively weak signature of emergence in (^Toi), when 
compared to the signature in {5rx cos 0) , and {Sry sin 0) , 
as this is not a simple converging fiow. Thus we may 
not be seeing a signal of emergence so much as a signal 
of supergranules. This in itself would be an interesting 
result, as it demonstrates a preference for emergence to 
occur at the boundary between supergranules. 

We note, however, the typical fiow speed for a super- 
granule of several hundred meters per second (Rieut ord] 
Rincon|2Q lQ) is an order of magnitude larger than our 
ensemble average of observed travel time shifts. This 
could be a result of imperfect alignment between the 
boundary and the emergence location, keeping in mind 
that our averaging disk is large compared to a typical 
supergranule size. By averaging over multiple supergran- 
ules for each emergence, the average travel time differ- 
ence would be greatly reduced. The difference in fiow 
speeds could also result from weakening fiow toward the 
edge of the sup ergranule. 

In [Paper III] we attempt to disentangle the relation- 
ships among some of the parameters where a difference 
between PE and NE samples is found. For example, do 
the Srmn maps contain in formation t hat is not in the 
maps of l- B^I? In addition, Paper III uses Discriminant 



Analysis (Kendall et al. 1983^ to determine which mea- 



surements arFlDesTllHetomstin the PE and NE 
populations, as a complement to the ensemble averaging 
performed here. 
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Second, when a subset of eleven "best" active re- 
gion candidates (emerging with a fast, monotonicahy- 
increasing flux history into an essentially field-free area) 
were considered, some of the differences described above 
disappeared. Specifically, there were no longer detectable 
j) or (Sroi). However, the direction- 
parameters continue 



differences in {Sr^ 



weighted {Stx cos 0) and {Sry sin i 

to show a statistically significant difference (of up to 2 s), 
for moderate depths. This subset includes the much- 



studied NOAA AR 10488 (see Paper I| for the full list). 

The results presented here place strong constraints 
on models of the emergence of AR. We have shown 
that in the ~ 24 hours before emergence that the sub- 
surface fiows, averaged over 6 hours, are no more than 
^ 15 ms~^. This raises the question of the time evolu- 
tion of the stron g (100 ms~ ^) retrograde fiows suggested 
by models (e.g. Fan 2008[ ): why have we not observed 
these fiows? how do these retrograde fiows interact with 
the near-surface shear layer? 

It is not clear how to reconcile the results presented 
here with the strong pre-emergence redu ction (of order 
10 s) i n the mean tr avel-time seen by Ilonidis et al. 



(2011). Ilonidis et al. ( ,2011, ) suggest that this signature 



may be caused by a rising tlux concentration crossing up- 
wards through the depth of 60 Mm, where their analysis 
is sensitive. Here, we use measurements that are sen- 
sitive to depths shallower than about 20 Mm, and see 
no signature in Srmn larger than a fraction of a second. 
Further measurements are needed to connect these two 
conclusions; why does the helioseismic impact of the ris- 
ing tube all but disappear as the tube approaches the 
photosphere? 

We have found that, on average for our pre-emergence 
sample, there are surface magnetic fields present at the 
emergence site but which vary only slowly (see Fig. 10 
of Paper I) over the day before fiux emergence begins 
in earnest. As discussed above, this may be due to a 
preference for emergence to occur at the boundaries be- 
tween supergranules where "quiet" magnetic fiux typi- 
cally accumulates continuously. As another possibility, 
we speculate that this may perhaps be a result of the 
interaction between convection and the rising magnetic 
fields, i.e., that the portion of the fiux that is caught in 
the fastest upfiows arrives at the surface well before the 
bulk of the fiux tube. Yet another possibility may be 
that fiux emergence occurs preferentially into remnant 
field. A more detailed analysis is needed to disentangle 
these effects. 

The data analysis we have presented here suggests that 
the r apid emergence process simulated by Cheung et al. 
([2010 ) is not typical. Horizontal fiows of order km s~^ 
extending over tens of Mm would produce signals in Stx 
and dry of order tens of seconds - well above our noise 
level. Notice, however, if fiows of this strength were to 
develop only after the emergence time, they would not 
appear in the current analysis. The same is true for the 
retrograde fi ows s e en in the rising fiux tube simulations 
discussed in Fan (2008). The picture we find h ere is 
apparently more compatible with the scenario of [Stein 



et al. (2011 ), in which weak field is brought to the surface 



by convection, which is itself only weakly altered by the 
magnetic field. 
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